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Abstract 


We study the stau lifetime in a scenario with the LSP taken to be a neutralino and the NLSP 
being a stau, based on the minimal supersymmetric Standard Model. The mass difference between 
the LSP and NLSP, 5m, must satisfy 5m/m ^ ~ a few % or less for coannihilation to occur, where 
m^ is the neutralino mass. We calculate the stau lifetime from the decay modes f —» yr, yz/ r 7T, 
and yz/ r A i{e)v^ and discuss its dependence on various parameters. We find that the lifetime is 
in the range ICY 22 TO 16 sec for 10 -2 < 5m < 10 GeV. We also discuss the connection with lepton 
flavor violation if there is mixing between sleptons. 
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I. INTRODUCTION 


The existence of non-barionic dark matter is now confirmed and its density has been 
quantitatively estimated Hi However its identity is still unknown. One of the most 
prominent candidates is the weakly interacting massive particle (WIMP) 0 0- El EJ- 

As is well known, the supersymmetric extension of the Standard Model provides a stable 
exotic particle, the lightest supersymmetric particle (LSP), if R parity is conserved. Among 
LSP candidates, the neutralino LSP is the most suitable for non-barionic dark matter since 
its nature fits that of the WIMP 0,0. Neutralinos are a linear combination of the supersym¬ 
metric partners of the U(l) and SU(2) gauge bosons (bino and wino) and the Higgs bosons 
(Higgsino). They have mass in a range from 100 GeV to several TeV and are electrically 
neutral. The lightest neutralino is stable if R parity is conserved. 

Since the supersymmetric extension of the Standard Model is the most attractive theory, 
the nature of neutralino dark matter has been studied extensively 0. In many scenarios 
of the supersymmetric model, the LSP neutralino consists mainly of the bino, the so-called 
bino-like neutralino. In this case, naive calculations show that the predicted density in the 
current universe is too high and it is necessary to find a way to reduce it. One mechanism 
to suppress the density is coannihilation ]10]. If the next lightest supersymmetric particle 
(NLSP) is nearly degenerate in mass with the LSP, the interaction of the LSP with the 
NLSP is important in calculating the LSP annihilation process. For coannihilation to occur 
tight degeneracy is necessary, since without coannihilation the LSP decouples from the 
thermal bath at T ~ m/20 |uj, where m is the LSP mass. Therefore the mass difference 
dm must satisfy 8m/m < a few %, otherwise the NLSP decouples before coannihilation 
becomes dominant. Furthermore, if the degeneracy is much tighter, we would observe a line 
spectrum of photons from pair annihilation of dark matter |l2|,[lj3|, since the annihilation 
cross section of dark matter would be strongly enhanced due to the threshold correction. 

A candidate for the NLSP is the stau or stop in many class of MSSM, and in this 
paper we study the lifetime of the stau-like slepton having mass degenerate with the LSP 
neutralino. For the neutralino LSP to be dark matter, very tight degeneracy is required in 
mass between the NLSP and the LSP neutralino. In particular the heavier the LSP is, the 
tighter the degeneracy must be 0]. For such a degeneracy, the NLSP is expected to have a 
long lifetime due to phase space suppression 14 Il5j|. 

An alternative scenario for a long-lived scalar particle is the gravitino LSP. Considerable 
work has been devoted to the long life NLSP in the context of the gravitino LSP. In this case, 
due to the small coupling between a superWIMP (including the gravitino) and the NLSP, 
the lifetime of the NLSP becomes very long LLj], 1 SB . To determine the most likely 
candidate for the LSP, we can accumulate and identify [UJ, |2jQ| the candidate for long-lived 
NLSPs and compare the nature of the particles including couplings. 

An implication on the fundamental feature of quantum mechanics would be brought by 
a long life stau. As noted in ref. [21], a Small-Q-value S-wave (SQS) decay can exhibit non¬ 
exponential decay. A small Q value naively implies a small mass difference. For a gravitino 
LSP, the decay must occur in a P wave since the gravitino has spin 3/2. In contrast, 
in our case the decay can be S wave since all daughter particles have spin 0 or spin 1/2. 
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Therefore, to extract fundamental parameters from observations we must take special care in 
interpreting the results. This study offers the opportunity to examine a fundamental problem 
of quantum physics in collider physics. Hence it is worthwhile studying the proposition that 
the stau is the NLSP and the lightest neutralino is the LSP, and that their masses are tightly 
degenerate. 

In Sec. ED we present the relevant Lagrangian and calculate the decay rate of the stau. 
In Sec. E3U we study the parameter dependences of the lifetime. We then consider the 
connection with lepton flavor violation (LFV) in Sec. ITVI Finally we summarize our results 
in Sec. El 


II. DECAY RATE 

In this section, we calculate the decay rate of the stau NLSP. The stau is a mass eigenstate 
consisting of superpartneres of left- and right-handed taus, 

f = cos 0 t ?l + sin 6 T e~ llT TR. (1) 

Here, 6 T is the mixing angle between 77 , and tr, and y T is the CP violating phase. The decay 
mode is governed by the mass difference, dm = ttinlsp—ttilsp; according to kinematics. That 
is, the lifetime of the stau depends strongly on dm. 

We consider the following four decay modes, f —> yr, t —> yny7r, f —> y/ikyz^, and 
r -> xev T v e (see Fig. Note that the NLSP can decay into other particles, for example, 
if dm > 1.86 GeV, a D meson can be produced in the stau decay but f —1 yr is dominant 
in this dm region since the D meson production process is suppressed by couplings and 
propagators. In the 3-body and 4-body decay processes, f —»■ yi/ r 7r, f —> xpu T u fll and 
f —»■ xeu T u e , diagrams can be formulated with charginos as intermediate states, however, 
such processes are strongly suppressed by the chargino propagator and we can safely ignore 
them Q. 



In this paper, we consider the small mass difference case and hence we can ignore the 
momentum in the W boson propagator. Thus the interaction Lagrangian is given by 


G AG 

Cint = t*x(9lPl + 9rPr)t + -^=u t x^PltJ^ + -^(/ 7 M PLZ/ / )(h T 7 At P L r) + h.c. 


( 2 ) 
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The first term describes stau decay into a neutralino and a tau. Here, Pl and Pr are the 
projection operators and 57 and g R are the coupling constants given by, for example in the 
bino-like neutralino case, 


9l = 7 =~~—~ sin 9 W cos 0 T , g R = -— sin 9 W sin d r e llT , (3) 

V 2 cos 6 *vi/ cos Ow 

where g is the weak coupling constant and 9w is the Weinberg angle. The second and third 
terms describe tau decays into a pion and/or leptons, where G is the Fermi constant. 

The detailed calculation of the stau decay rate is given in the appendix El Here we show 
approximate formulae. In the region dm > m T , 2-body decay (see Fig. Ha)) is allowed 
kinematically and is dominant. The decay rate of the 2-body final state is approximately 

r 2 -bod y = ~ m r (( 9l + \9r\ 2 )& m ~ 2 Re[g L g R }m T ) , (4) 

where m f: m^, and m T are the masses of f, x, and r, respectively. 

The 2-body decay is forbidden kinematically for dm < m T , and pion production (see 
Fig- Hb)) is dominant if dm is larger than the pion mass m,. The pion production 3-body 
decay rate has the approximate form 


3-body 


G 2 f 2 cos 2 9 C 


((dm) 2 -ml) 


5/2 


X 


210(27r) 3 m^m^ 

g\dm (4 (dm) 2 + 3 m 2 ) — 2 Re[gLg R }m T (4 (dm) 2 + 3 ml) + 7\g R \ 2 m 2 dm 


( 5 ) 


Here f n is the pion decay constant and 9 C is the Cabbibo angle. 

Incidentally, we note that a quark cannot appear alone in any physical processes, a point 
that was missed in ref. [14j. Hence, u and d quarks appear only as mesons and the u, d 
production process is relevant only for dm > m n . 

Finally, when the mass difference is smaller than the pion mass, 4-body decay processes, 
f —> xg,v T v^ and f —> x ez/ r^ e , are significant (see Fig. He)). The approximate decay rate is 
calculated as 


G 2 


((dm) 2 — m 2 ) 5/2 


4-body 945 ( 27 r) 5 m^mjl 

2 gl(dm) 3 (2 (dm) 2 — 19m/) — ARe{gLg R \m T (dm) 2 (2 (dm) 2 — 19m/) 
+ ?>\g R \ 2 m 2 dm (2(dm) 2 — 23m/) 


x 


( 6 ) 


Here mi is the charged lepton (e or /i) mass. 


III. PARAMETER DEPENDENCE 

In this section, we discuss the parameter dependence of the stau lifetime and the cosmo¬ 
logical constraints for the parameters in the bino-like LSP case. From Eqs. (ED, ©, and ©, 
we see that the stau lifetime depends on dm,9 T ,m and 7 T . 
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A. m ^ 

First, we examine the neutralino mass dependence of the lifetime. We consider the 
cosmological constraints for dark matter to get the mass range of the LSP. It is well known 
that the dark matter relic density is reduced by the coannihilation process and hence the 
neutralino mass can be heavier than it would be without coannihilation. Accounting for the 
coannihilation process gives a neutralino mass range (the first ref. of [f|) of 

200 GeV < m* < 600 GeV. (7) 

Here, we use CMSSM bound as a reference value of though we study the stau lifetime 
in general MSSM framework, since it is not strongly dependent on m^ as noted below. This 
is consistent with the cosmological constraint, 0.094 < < 0.129. 

It is clear from Eqs. and & that the stau lifetime is proportional to the neutralino 

mass and therefore in the mass range given in Eq. (ED, the stau lifetime varies by a factor 
of three. Taking this into account, we use only = 300GeV in our figures. 


B. dm 


Next, we consider the dm dependence of the stair lifetime. For coannihilation to occur, 
the mass difference must satisfy dm/m ^ ~ a few % or smaller lo||. Therefore, dm must 
be smaller than a few GeV. The dm dependence of the total stau lifetime and the partial 
lifetimes of each decay mode in this dm region is shown in Fig. [21 In Fig. [21 we set the values 
as follows: m * = 300 GeV, 9 r = 7t/ 3, and y T = 0. 

^Frorn Fig. [21 we can see which mode is dominant for a certain dm. The figure shows that 
the lifetime increases drastically as dm becomes smaller than the tau mass. This is because 
taus are produced in the region dm > m T , while in the region dm < m T , taus cannot be 
produced and instead pions appear in the final state. At dm = m n the lifetime increases 
slightly. This is due to the fact that the dominant mode changes from 3-body to 4-body 
decay. In contrast, at dm = m^ the lifetime does not increase much, even though above 
this mass, muons can be created. This is because at the pion mass, the muon production 
process is already kinematically suppressed and the electron production process governs the 
stau decay. 

To understand the dm dependence of the lifetime quantitatively and intuitively we de¬ 
termine the power of dm in the decay rate, considering stau decay into neutralino and n — 1 
massless particles. The dm dep endence of the decay rate is determined by the phase space 
and the squared amplitude Jl4j. 

First, we examine the dm dependence by considering the phase space. For 2-body decay, 
a phase space consideration gives 


d<jF> = 


dQ 

32vr 2 


1 - 


m x 


m £ + dm , 


oc dm . 


( 8 ) 
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FIG. 2: Total lifetime and partial lifetimes of each decay mode as a function of dm. The lines label 
electron, muon, pion, and tau correspond to the processes f —» yezy-zy, f —► xnv T v tl , t —» x^tT, 
and f —> xt, respectively. Here we take m% = 300 GeV, 0 T = ir/3, and y T = 0. 


By using a recursion relation between dq^ and d(j)^ n the phase space of n-body decay 
renders the Sm dependence as 

/ 5m 

oc (Sm) 2 ^ n ~ 2 ^ +1 . (9) 

Second, we consider the dm dependence from the squared amplitude. If all of the n — 1 
massless particles are fermions, the squared amplitude depends on 5m as 

M (n) oc (dm) 71 - 1 , (10) 

since it depends linearly on the massless fermion momentum. Thus, we obtain the depen¬ 
dence of the decay rate on dm for a final state of only fermions, 

r W oc M (n) x dc/> (n) oc (dm) 3n ~ 4 . (11) 

In contrast, if one pion(NG-boson) appears in the stau decay process, the dm dependence 
of the squared amplitude becomes 

M (n) oc ( dm) n . (12) 


6 




























This change in the decay process is due to the fact that the amplitude of the pion production 
is proportional to the pion momentum. Namely, the squared amplitude of the pion produc¬ 
tion process is proportional to the pion momentum squared. Thus the 8m dependence of 
the process, in which one pion is involved, is 

T (n) oc (8m)^ n ~ 1] . (13) 

In the massless limit of external line particles, the 8m dependences of our results, calculated 
in appendix El are 

T2-body OC {8Til') , 

r 3 -body OC ( 8m ) 6 , 

r 4 -body OC (8m) 8 , (14) 

which are consistent with Eq. (ED and Eq. ED- 

More precisely, by taking into account the masses of the produced particles, we get the 
complicated 8m dependences 

IYbody OC (8m) {(8m) 2 - m 2 r ) 1/2 , 
r 3 -body OC (8m) {(8m) 2 - m 2 ) 5/2 , 

r 4 - bod y oc (8m) 3 {(8m) 2 - m 2 ^' 2 . (15) 

These equations clearly show the 8m dependence of the stau lifetime. 

C. 7t 

We next consider the 7t dependence of the stau lifetime. Figure 01 shows the lifetime as 
a function of the CP violating phase. We set the other parameters to be m ^ = 300GeV, 
8m = 0.5GeV, and 0 T = tt/3. 

^Frorn Fig. 01 it is clear that the CP violating phase does not greatly affect the stau 
lifetime, and so we fix 7t = 0 (no CP violation). 

As expressed in Eq. (IbD . the effect of CP violation appears in the Re[gLgn\ terms only. 
Since the coefficients of the Re[gLgn\ terms are smaller than those of Ig^J 2 , it is again clear 
that the CP violating phase does not greatly affect the stau lifetime. 


D. e T 

The 9 t dependence of the stau lifetime is as strong as the 8m dependence. Figure 0] shows 
the 6 r dependence of the stau lifetime with m% = 300GeV and 7r = 0. 

We can see from the right panel in Fig. 01 that for 8m -C m T , tr decays much more 
quickly than ?l . This can be understood by considering two steps. First, we note that 
only left-handed virtual taus contribute to the final state u T . Second, Tr converts to tl by 
picking up m T in the tau propagator, while tr converts by picking up the momentum p T in 
the propagator. Since p T ~ 8m -C m T , the former contribution is much larger and hence 
there is a strong dependence on 0 T . 
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FIG. 4: Stau mixing angle and stau lifetime with m ^ = 300GeV and y T = 0. Left panel: 5m = 2.0 
GeV, right panel: 5m = 0.2. 


E. Constraint from Big Bang Nucleosynthesis 


Finally, we consider the constraints from Big Bang nucleosynthesis (BBN). The stau 
lifetime will be constrained by the standard BBN scenario. Particularly, in the case where 
a stau decays into hadrons, the constraint on the stau lifetime is more stringent j22[ Esj. 
At t > 100 (sec) (t is the age of the universe), pions emitted from a stau decay also decay 
before they interact with the background nucleons, so they do not affect BBN. However, if 
pions are emitted at t < 100 (sec), they inter-convert a background proton and neutron into 
























each other, even after the normal freeze-out time of the n/p ratio. Therefore we must take 
careful consideration to the parameter region in which pions are emitted for (BBN start) 
<t< (BBN end), i.e. 10“ 2 (sec) < t < 10 2 (sec). Naively interpreting Fig. 38 in ref. jQ, 
we can conclude that the corresponding parameter region can neglect the constraints from 
BBN. This is inferred from the fact that the f number density is roughly two orders lower 
than that of baryons, Y f ~ 10 —12 , since rriy, is of the order of several hundreds of GeV, 
and the released energy is less than 1 GeV. ffere Y f is the ratio of the stau number density 
to the entropy density. Thus E vis Yf < 1CT 12 GeV, which does not influence the BBN for 
t < 10 2 sec. Incidentally, it is apparent from Fig.5 that stau lifetime can not exceed the age 
of universe (about 4.1 xlO 1 ' sec). It is consistent with the fact that an exotic heavy charged 
particle has not been observed yet. 


[X7t X 



5m (GeV) 


FIG. 5: Contour plot of stau lifetime in sec as a function of 5m and d T . To emphasize the BBN 
era, we draw thick lines for 10~ 2 and 10 2 sec. 


IV. CONNECTION WITH LEPTON FLAVOR VIOLATION 

We next consider LFvj29|. The NLSP slepton might be a linear combination of flavor 
eigenstates. It is expected that lepton flavor violating events will be observed due to this 
mixing, such as /i —> ey. If we observe r —> ey or r —> /iy events, then within the context of 
the minimal supersymmetric Standard Model (MSSM), we would conclude that the selectron 
forms part of the NLSP: 

0NLSP = Nie + \J 1 — lV, 2 f. (16) 

The branching ratio of r —> ey is roughly proportional to Nf. The current upper bound on 
the branching ratio, < O(10 -7 ), gives a poor constraint, at most A’i < 0.1. If N 1 ^ 0, the 
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NLSP slepton can decay, 


0NLSP —> X + G 


( 17 ) 


with a lifetime of 


1.666 x 10" 18 


m x \ 
300GeV/ 


/ 1.0GeV\ 2 
\ 5m ) 



(18) 


We may interpret the result as implying that the NLSP slepton consists only of a scalar 
electron since it decays only into an electron. The lifetime is so short, even if 5m < m T , that 
we cannot accumulate the NLSP. However if we can measure the lifetime, we can determine 
the fraction of the scalar electron in the NLSP slepton 
violating processes of both charged leptons 


m 


3pton. Furtl 

3, 0 0 


Further evidence is given by flavor 
and the neutrino |28l |. 


It is 


important to interpret all the LFV processes together. 


V. SUMMARY 

We have studied an MSSM scenario in which the LSP and NLSP are a bino-like neu- 
tralino and a stau, respectively. Since the mass difference is, in many cases, assumed to be 
degenerate, from the requirement of coannihilation, we paid special attention to the very 
small 5m case. We calculated the partial lifetimes for the decay modes shown in Fig. 0] 

We have investigated the stau lifetime dependence on 5m, d T , 7 T , and m^, considering 
cosmological constraints. The lifetime strongly depends on 5m and 0 T , while it is almost 
independent of y r and m*. 

The stau lifetime dependence on 5m changes as each threshold is crossed. When 5m 
is larger than m T , the lifetime increases in proportion to (<5m)~ 2 as 5m decreases. In the 
range m T > 5m > m n the lifetime obeys the scaling ~ (<5m) -6 . Below it grows with 
(5m)~ 8 . The 5m dependence of the stau lifetime can be largely understood by counting the 
mass dimension of phase space and the squared amplitude in the massless limit of Standard 
Model particles. While the massless limit is a good approximation in regions far from the 
thresholds, the 5m dependence near the thresholds are given by Eq. (USD- 

The lifetime also strongly depends on 9 T , as shown in Fig. |3 tr contributes to 3- and 
4-body decay processes by picking up the m T term in the intermediate r propagator, while 
tl picks up the p T term. Since p T ~ 5m <C m T , the contribution for tr is much larger and 
hence there is a strong dependence on 9 T . 

As is seen in Fig. |21 if 5m is smaller than m T , a stau can be very long-lived. This fact 
may need to be taken into account in studies at the Large Hadron Collider (LHC) and 
the International Linear Collider (ILC). The super WIMP case is similar. However, in the 
gravitino scenario an energetic r is apparently produced, while in our scenario we should 
observe very low energy n, /i, or e. Therefore a clean experiment is required. The ILC would 
be most suitable for investigating the nature of the NLSP slepton. 

Another candidate for the NLSP is the stop. However, it is more complicated to investi¬ 
gate this possibility. A stop must always be accompanied by quark(s); it forms a meson-like 
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fermion with one quark or a baryon-like boson with two quarks. It is almost impossible to 
calculate the exact mass eigenvalue using QCD. It is beyond the scope of this paper and is 
left for future work. 

We have also discussed lepton flavor violation due to slepton mixing. If there is even a 
tiny component of a scalar electron or a scalar muon in the NLSP “stau”, the decay signal of 
the NLSP will be completely different from the pure stau case. The NLSP slepton undergoes 
2-body decay into the accompanying electron or muon. Since it is a 2-body process, it occurs 
very quickly, ~ (10 _20 )iVf 2 sec where N\ represents the portion of the scalar electron or 
scalar muon, as shown in Eq. USD- As this mixing causes charged/neutral lepton flavor 
violation, it is very important to compare the NLSP slepton decay with other processes such 
as r —*■ e(/i) 7 - To fully clarify the nature of the NLSP, we need to interpret all the LFV 
processes together. 

SQS decay takes place in a case of small mass difference. It will appear strongly, if decay 
rate is small. Since long lifetime means small decay rate, long-lived stau will give a good 
chance to observe a non-exponential decay. In contrast, if there is no LFV in slepton mixing, 
the pure stau has a very long lifetime and it is then possible to experimentally observe SQS 
decay [21]. This is important to check the fundamental feature of quantum mechanics. 

Hence, the small 5m case is very interesting and its study is very important. 
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APPENDIX A 

In this appendix, we show the exact and approximate stau decay rates. We use the exact 
decay rate in formulating the figures. We make the following approximations: we keep only 
the leading order term of (5m/m^) and we replace the denominator of the r propagator by 
m 2 in the 3- and 4-body cases. We calculate the decay rates of the three processes shown 
in Fig. |T| with the Lagrangian of Eq. (11- 


2-body decay 

The decay rate of the 2-body decay process (see Fig. H*a)) is given by 
1 


r2-body = 16 3 (ra* + m\ + m* - 2m~m~ - 2m>; - 2m\m/)* 


x {(3l + I- "4 - rnl) - 4:Re[g L g R ]m T m *} . (Al) 
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For the analysis discussed in Sec. El we approximate the decay rate as 

1 


2-body 


Anm^ 


^ {5m) 2 - m 2 T (( g\ + \g R \ 2 )Sm - 2Re[g L g R ]m T ) . 


(A2) 


3-body decay 


The decay rate of the 3-body decay process (see Fig. m) is calculated as 
G 2 ff cos 2 e c ((dm) 2 - ml) 


3-body 


647r 3 ml 


x J dx^[(5m) 2 — g 2 ) ((5m + 2m^) 2 — q 2 ) — 


,2 _ m 2y _|_ ( mr r T y 


x (?/ - m l) 


t (dijQf + \g R \m T )((5m) + 2m^5m - g f ) - Re[g L g R ]m^m T q f 


(A3) 


Here q 2 is given as 

q 2 = (5m) 2 — [(5m) 2 — m 2 ) x , (A4) 

where the index /(= 7r, e, /x) denotes a massive particle, except the neutralino, in the final 
states; f — 7r in the 3-body case. T t is the tau decay width and (m T r r ) 2 is added to the 
denominator of the tau propagator for the region 5m > m T . 

The approximate decay rate is 


r 


3-body 


G 2 fl COS 2 0 C 
210(27t ) 3 m^ml 


[(5m) 


m: 


5/2 


X 


g\5m (4 (5m) 2 + 3m 2 ) - 2Re[g L g R }m T (4 (5m) 2 + 3m 2 ) + 7| g R \ 2 m 2 T 5m 


(AS) 


4-body decay 


In the 4-body decay processes (see Fig. me)), the decay rate is given by 
G 2 ((5m) 2 — mf) 


4-body 


24(27r) 5 m| 

X l dx^({SmY- q f, ((5m + 2mtY- q )) (nav) 

1 


2 n A 

q f 


X 


4 (0l 9/ + I Shilov) ((5m) + 2my,5m - q f ) - Re[g L g R ]m^m T q f 
x <j 12mfgJ log + ( g J - mf)(q 4 f - 8m z 2 gj + mf) )> , (A6) 


mf 


where l — e, n and qf is given by Eq 
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We can approximate the decay rate as 

= 945 ( 2 ^, ((6m ? ~ m V' 2 

x 2g 2 (hm) 3 (2 {5m) 2 — 19m 2 ) — ARe[g L g R ]m T (8m ) 2 (2(<5m) 2 — 19m 2 ) 
+ 3\g R \ 2 rn 2 5m (2(hm) 2 — 23m 2 ) 
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